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ABSTRACT 

Most late-type stars in the solar neighborhood have velocities similar to the local standard 
of rest (LSR), but there is a clearly separated secondary component corresponding to a 
slower rotation and a mean outward motion. Detailed simulations of the response of a 
stellar disk to a central bar show that such a bi-modality is expected from outer-Lindblad 
resonant scattering. When constraining the run of the rotation curve by the proper motion 
of SgrA* and the terminal gas velocities, the value observed for the rotation velocity 
separating the two components results in a value of (53 ± 3) kms" 1 kpc^ 1 for the pattern 
speed of the bar, only weakly dependent on the precise values for R and bar angle <fi. 

Subject headings: Galaxy: kinematics and dynamics - - Galaxy: structure - - solar 
neighborhood 
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1. Introduction 

Even though wc now know beyond reasonable doubt that 
the Milky Way is barred, the structure of the bar, its ori- 
entation and pattern speed are still subject to substantial 
debate, mainly because of the edge-on view and dust ob- 
scuration. Consequently, the properties of the bar have 
been inferred rather indirectly from IR photometry ( Blit2 
k Spergel 199lt |Dwek et al. 1995| ; |Binney et al. 1997|) , 
asymmetries in the distribution or magnitud e of stars 
flWhitelock fc Catchpole 1992j; [Weinberg 1992|; |Nikolaev| 
L Weinberg 19971 ; |Stanek 1995| ; |Sevenster 1996|; [Btanek 



3t al. 1997), and gas velocities (Liszt fc Burton 1980 



Binney et al. 1991 ; Englmaier fc Gerhard 199S ; Weinei 
fc Sellwood 1999 ). While all these studies are based on 
observations of the inner Galaxy itself, it is the goal of 
this letter to constrain the properties of the bar by its 
effect on the stellar velocity distribution observed in the 
solar neighborhood. 

2. The Local Stellar Velocity Distribution 

In contrast to the structure of the inner Galaxy, the lo- 
cal stellar velocity distribution can be observed to much 
greater detail than for any external galaxy. Figure ^ 
shows the distribution /(it, v) in radial (u) and rotational 
(v) velocities (relative to the LSR) that has been inferred 
from the Hipparcos data ( |ESA 1997| ) of 6 018 late-type 
stars dDehnen 1998| fc 1999a). 

This distribution shows a lot of structure, most of 
which is real (noise has been suppressed). Apart from 
an overall smooth background and a low-velocity re- 
gion (solid ellipse in Fig. [I]), within which / has several 
peaks associated with moving groups, there is also an 
intermediate-velocity structure (broken ellipse) contain- 
ing ~ 15% of the late- type stars but hardly any early- 
type stars, which mainly populate the moving groups 



(Dehnen 1998; Chereul et al. 1998; Asian et al. 1999 



Skuljan et al. 1999| ). This secondary component, which 
consists predominantly of stars moving outwards (u < 0) 
and is clearly separated from the low-velocity peak, was 
only vaguely recognizable in pre-Hipparcos data sets; 
the resulting mean outward motion of stars with v <J 
— 30kms _1 is also known as u-anomaly. 




u [km s" 1 ] 

Fig. 1. — The distribution f(u,v) inferred from Hipparcos data 
for late-type stars (3 527 main-sequence stars with B — V > 0.6 
and 2 491 mainly la te-type non-rn airi-sequence stars, high- velocity 
stars excluded, see [Dehnen 199£ ). u and v denote the velocities 
towards I = 0° and £ = 90° w.r.t. the LSR measured by Dehnen 
& Binney (1998); indicates the solar velocity. Samples of early- 
type stars contribute almost exclusively to the low-velocity region 
(solid ellipse). The region of intermediate velocities (broken ellipse) 
is mainly represented by late-type stars, ~ 15% of which fall into 
this region. 

the intermediate- velocity stars and for their clear separa- 
tion from the low-velocity stars. The first fact rules out 
any axisymmctric equilibrium, while the second strongly 
hints towards an orbital resonance. 

3.1. The Bar's Outer Lindblad Resonance 

Such a resonance is provided by the Galactic bar: the 
outer Lindblad resonance (OLR), which occurs when 



^0 



(1) 



Here, and ujr are the azimuthal and radial orbital fre- 
quencies, while fib is the pattern speed of the bar. After 
each radial lobe, a star on a resonant orbit is at the same 
position w.r.t. bar, and hence, is always pushed in the 
same direction forcing it onto a different orbit. According 
to the modeling of the gas motions in the inner Galaxy 
([Englmaier fc Gerhard 1999|; |Weiner fc Sellwood 1999|), 



3. What Causes this Bi-Modality? 



the radius -R lr corresponding to the OLR of circular 
orbits is not far from Rq. 

Orbits inside (outside) the OLR, i.e. with uj^-\- \ur > 
fib (< f^b), are elongated perpendicular (parallel) to the 
bar's major axis and are moving on average outwards 
Attempts to explain this seeming anomaly by non-equilibriun^ inwards) for bar angles <p G [0°,90°]. The bar angle is 
dynamics, like the dispersion of a stellar cluster or the defined as the azimuth of the Sun w.r.t. the bar's major 
merger of a satellite, fail because of the exclusively late axis and, according to IR photometry, lies in the range 
stellar type, disk- like kinematics (v nearly as LSR), 15° to 45°. 

and rather high metallicity ( Raboud et al. 1998| ) of the Thus, the OLR of the Galactic bar can naturally ex- 



intcrmediate-velocity stars. Thus, dynamical equilib- 
rium must account both for the mean outward motion of 



plain the low- versus intermediate-velocity bi-modality 
of f(u, v), provided the Sun is outside the OLR ( Dehnen 
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1999a). In this picture, the depression between the 
modes of f(u,v) is due to exactly resonant orbits, and 
from the observed w-velocity of 



= (-31 ±3)kms" 



(2) 



for the saddle point between the two modes, we can es- 
timate that Rq — -Rolr — \vo-lb./vo\Ro — lkpc (for a 
fiat rotation curve). Thus, only stars with epicycle am- 
plitudes ;> lkpc, i.e. predominantly late-type stars, can 
visit us from inside the OLR, in agreement with the ab- 
sence of early- type stars in the secondary component. 

This proposal is in line with the ideas of Raboud et al. 
(1998), who observed a high metallicity and outward 
mean motion for disk stars with v < — 30kms _1 and 
interpreted this as evidence for the Galactic bar, since 
in Fux's (1997) iV-body simulations of the Milky Way 
u < at the solar position. However, these authors did 
not relate the effect to the OLR nor were they able, for 
lack of resolution both in the data and the model, to 
make detailed quantitative investigations. 

3.2. The Response of the Disk to a Stirring Bar 

In order to quantify the behavior of a warm stellar disk in 
presence of a central bar, I performed numerical simula- 
tions in which a bar is grown in an initially axisymmetric 
exponential disk. These simulations, presented in detail 
in Dehnen (1999b), use a technique, which is similar to 
that of Vauterin & Dejonghe (1997) and allows for very 
high resolution. 

The simulations show that a second mode appears in 
f{u, v) at many positions in the Galactic disk, provided 
the velocity dispersion is high enough. The velocity v OLR 
of the division line between the modes depends essen- 
tially on four parameters: (1) our distance to the OLR, 
quantified by the ratio R lr/Ro', (2) the bar angle (/>; the 
(3) shape and (4) normalization of the rotation curve, 
quantified by (3 = dlni> c /dlni? and vq = v c (Rq). The 
strength of the bar, its morphology, and the details of 
the stellar DF hardly affect the velocity v OLR , but may 
change the strength of the secondary component. 

The simulations used rotation curves v c cx R and a 
simple quadrupole to model the non-axisymmetric con- 
tribution of the bar potential. Figure || shows the result- 
ing f(u, v) for /3 = 0, <f> = 25°, and three different values 
for f2b, which is related to -R lr by 



fioLR _ / ^0 

~rT ~ l Ob 



1/(1-/3) 



(3) 



with f^o = vo/Rq- The secondary component is weaker 
and occurs at larger \v\ the farther the OLR is from the 




Fig. 2. — Simulated f(u,v) after the growth of a central bar. 
The initial axisymmetric disk has scale length 0.33i?o an d velocity 
dispersion reminiscent of the old stellar disk. The rotation curve 
is assumed fiat and the bar angle is fixed at <f> = 25° . The pattern 
speed decreases from top to bottom panel. 

Sun. For 15° < < 45°, -0.2 < (3 < 0.2 and 0.8 < 
Rolr/Rq < 0.95, i>olr is well approximated by 



"OLR 



1 + /3 
'1^8 



t'o 



Ro Ot 



(& + e/3)v ,(4) 



+ v /(l + /3)/2_ 
where the coefficients a, b, and c are given in Table |. 

4. Implications for the Inner Galaxy 

We now assume that the observed low- versus intermediate- 
velocity bi-modality of /(it, v) is caused by scattering off 
the OLR of the Galactic bar. Then the observed velocity 
w OL r (||) constitutes a new constraint on the structure of 

Table 1: Best-fit values for (a,b,c) in equation (Q) 





a 


b 


c 


15° 


1.3549 


0.0761 


0.1362 


20° 


1.2686 


0.0642 


0.1120 


25° 


1.2003 


0.0526 


0.0892 


30° 


1.1424 


0.0406 


0.0711 


35° 


1.0895 


0.0298 


0.0538 


40° 


1.0420 


0.0200 


0.0423 


45° 


1.0012 


0.0103 


0.0316 
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Table 2: 


Adding the new constraint: 


fitting Vteim for \l 


> 30°, /* 


of Sgr A*, and v 


OLR • 


Rq - 


: 7, v = 198 ± 5.6, = -0.147 ± 0.025 


Rq = 


-- 8, v = 227 ± 6.4, p = -0.020 ± 0.022 


<t> 


fi b 


Rolr/Ro 


Rcb/Rq 


<t> 


fi b 




Rcr/ Rq 


15° 


51.4 ±1.4 


0.920±0.011 


0.594 


15° 


50.5 ±1.4 


0.954±0.010 


0.567 


20° 


52.1 ±1.5 


0.909 ±0.012 


0.586 


20° 


51.1 ±1.4 


0.943±0.010 


0.561 


25° 


52.8 ±1.5 


0.898 ±0.012 


0.579 


25° 


51.8 ±1.4 


0.932 ±0.010 


0.554 


30° 


53.7±1.5 


0.886 ±0.013 


0.572 


30° 


52.5 ±1.4 


0.920±0.010 


0.547 


35° 


54.5 ±1.6 


0.874±0.013 


0.564 


35° 


53.2 ±1.5 


0.908±0.011 


0.540 


40° 


55.3 ±1.6 


0.863 ±0.013 


0.557 


40° 


53.8 ±1.5 


0.897±0.011 


0.533 


45° 


56.2 ±1.6 


0.851 ±0.013 


0.549 


45° 


54.5 ±1.5 


0.885±0.011 


0.526 


Rq = 


= 7.5, v = 212 ± 6, = -0.080 ± 0.024 


Rq = 


= 8.5, v = 241 ± 6.8, P = 0.033 ± 0.021 


<t> 


fi b 


Rolb/Ro 


Rcr/Ro 


<t> 




Rolr/R() 


Rcr/ Ro 


15° 


50.8 ±1.4 


0.939 ±0.010 


0.581 


15° 


50.4 ±1.4 


0.967 ±0.009 


0.552 


20° 


51.5 ±1.4 


0.928±0.011 


0.575 


20° 


51.0 ±1.4 


0.956 ±0.009 


0.546 


25° 


52.2 ±1.4 


0.917±0.011 


0.568 


25° 


51.6 ±1.4 


0.944 ±0.009 


0.539 


30° 


52.9 ±1.5 


0.905±0.011 


0.560 


30° 


52.2 ±1.4 


0.932 ±0.010 


0.532 


35° 


53.7 ±1.5 


0.893 ±0.012 


0.553 


35° 


52.8 ±1.4 


0.921 ±0.010 


0.526 


40° 


54.4 ±1.5 


0.882 ±0.012 


0.546 


40° 


53.5 ±1.5 


0.909 ±0.010 


0.519 


45° 


55.2 ±1.6 


0.870±0.012 


0.539 


45° 


54.1 ±1.5 


0.898 ±0.010 


0.513 



the inner Galaxy, the implications of which we will now 
investigate. 

4.1. Other Constraints 

Unfortunately, the value of w lr does not depend on just 
one, but several parameters. Thus, in order to under- 
stand the implications of the observed value, we must 
combine the new constraint with other, independent con- 
straints. 

Determinations of the proper motion of Sgr A*, the 
radio source associated with the Galactic center, yield 
(in kms-^pc" 1 ) m = -28.0 ± 1.7 (|Rcid et al. 1999[ ) 
and fi e = -29.3 ± 0.9 ( packer fe Sramek 1999[ ), while 
the latitudinal proper motion \x\, is negligible. Sgr A* 
is generally thought to be associated with a black hole 
of mass - 2.6 x 1O 6 M ( |Eckart fe Genzel 1997| , |Ghcz 



;t al. 1998 ), which because of energy equipartition must 
be (nearly) at rest w.r.t. its surroundings. Thus, the 
observed proper motion is entirely due to the reflex of 
our own rotation around the Galaxy, and we obtain 



fi = (29.0±0.8)kms- 1 kpc~ 1 - v Q /R , 



(5) 



where u© = (5. 25 ± 0.62) km s 1 is the Sun's u-motion 



while non-axisymmetric effects are likely to add small 
deviations to this relation. 

4.2. Adding the New Constraint 

We now put the various constraints together by (i) as- 
suming (reasonable) values for Rq and (j), (ii) computing 
vq = SIqRq using (||), (iii) fitting P to the observed ter- 
minal velocity curve, and (iv) deriving the bar's pattern 
speed fib from the constraint on i> OL r whereby using the 
approximation (||). Note that the bar angle becomes im- 
portant only at stage (iv), i.e. vq and P are already fixed 
once Rq is chosen. 

For Rq = 7, 7.5, 8, and 8.5 kpc and values of (f> in the 
range 15° to 45°, which is required for the bar's mor- 
phology to be both similar to that of external galaxies 
and consistent with IR photometry, Table [2] lists the re- 
sults of this procedure when fitting wtorm at \£\ > 30° 
(the results hardly depend to this limit). 

An inspection of the table reveals that the derived 
value for fib is only weakly dependent on the assumed 
values for Rq and <f>: fib varies by about 10% with larger 
<fi corresponding to larger fib- Owing to this weak de- 
pendence, the pattern speed of the bar is rather tightly 
constrained even without accurate knowledge of (Rq, 4>): 



w.r.t. the LSR ( |Dehnen fc Binney 1998j ) 

In order to constrain the shape of the rotation curve, 
i.e. the parameter P, we use the terminal gas velocities 
inferred by Malhotra (1994, 1995) from various observa- 
tions in Hi and CO. For axisymmetric gas motions, the 
terminal velocity is related to the rotation curve by 



fib = (53 ± 3) km s 1 kpc 



(7) 



Table y also lists estimates for R ol , r /Rq obtained from 
equation (|^) and for i? CR /i?o obtained from the terminal 
velocities by the relation 



v c (Ro\sm£\) - v c (R )\sin£\, 



(G) 



(0 



sin^ 



(8) 



|sin£| = Rc R /flo 
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An error is not given, but the uncertainty is about 10%. 

Amongst the three constraints used to derive fib the 
biggest uncertainty comes from the proper motion of 
SgrA* (and its interpretation as reflex of the solar mo- 
tion). If in equation (|s|) one uses the value reported 
by Reid et al., instead of a weighted mean with Backer 
& Sramek's measurement, the resulting fib is about 2% 
smaller, while R C r/Ro hardly changes. 

The estimate (|7]) is also subject to systematic uncer- 
tainties, due mainly to the idealized model. Future sim- 
ulations which allow for a more general Galactic rotation 
curve and bar potential and are fitted simultaneously to 
the gas motion of the entire inner Galaxy and the lo- 
cal stellar velocity distribution are likely to reduce these 
uncertainties considerably. 

5. Conclusion 

The simulati ons presented in Dehnen (1999b) and briefly 
reported in §3.2 show that the bar's OLR inevitably cre- 
ates a bi-modality in f(u,v) over a wide range of posi- 
tions in the Galactic disc. The low- versus intermediate- 
velocity bi-modality in f(u,v) observed for late-type 
stars (Fig. |l|), has the same characteristics as these OLR 
features and is thus best identified as such. 

This in turn allows the observed v-velocity seperat- 
ing the two modes to be used as a new constraint on 
the parameters of the inner Galaxy, notably the rota- 



eds. Merritt, Va 



Sell- 



tion curve, the bar's pattern speed and projection angle. 



In conjunction with constraints for the rotation curve 
from the terminal gas velocities and the proper motion 
of SgrA*, this results in a value for the bar's pattern 
speed of (53 ± 3) kms -1 kpc -1 . The systematic uncer- 
tainty is likely similar in size to the internal error; both 
shall be reduced in more elaborate future modeling. 

The values for the corotation radius implied by this 
pattern speed and the terminal gas velocities are between 
0.5 and 0.6 times Rq, i.e. larger than the upper limit of 
0.5 given by Englmaier & Gerhard (1999), but consistent 
with Weiner & Sellwood (1999). Given that IR photome- 
try suggests a bar length i?b of about 0.45i?o, this implies 
a ratio R CR /Rb of about 1.25 ± 0.2, in agreement with 
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